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Recent experiments have introduced a new concept for 
analyzing the photoemission spectra of correlated electrons - 
the remnant Fermi surface (rFs), which can be measured even 
in systems which lack a conventional Fermi surface. Here, we 
analyze the rFs in a number of interacting electron models, 
and find that the results fall into two classes. For systems 
with pairing instabilities, the rFs is an accurate replica of the 
true Fermi surface. In the presence of nesting instabilities, 
the rFs is a map of the resulting superlattice Brillouin zone. 
The results suggest that the gap in Ca2Cu02Cl2 is of nesting 
origin. 

Recently, a new experimental tool has been introduced 
JIJ to parametrize photoemission (PE) data in strongly 
correlated metals: the 'remnant Fermi surface' (rFs). 
This is the locus of points in fc-space where the PE inten- 
sity associated with a particular quasiparticle peak falls 
to half of its peak value. For an ordinary metal, these 
points would correspond to the true Fermi surface, but in 
strongly correlated metals the points do not necessarily 
fall at the same energy - the rFs may display a consid- 
erable dispersion. 

Ronning, et al. [Q measured the rFs of Ca2Cu02Cl2 
(CCOC), a half filled Mott insulator, and compared 
it with the rFs of optimally doped Bi 2 Sr 2 CaCu208 
(BSCCO). When underdoped, the dispersion of BSCCO 
evolves toward that of CCOC 0. Qualitatively, the 
rFs of CCOC seems consistent with Luttinger's theorem, 
even though it displays a considerable dispersion. De- 
spite this, the rFs's of the two materials are strikingly 
different, and cannot evolve from each other via rigid 
band filling (they cross). A proper understanding of the 
rFs could lead to an improved model for the pseudogap 
in these materials. 

In this Letter we analyze the rFs expected for a vari- 
ety of interacting electron systems, and show that they 
do not necessarily provide information about the 'true' 
Fermi surface. The results fall into two classes, depend- 
ing on whether the interaction can be characterized as 
'nesting' or 'pairing'. Only in the latter case is the rFs 
a reliable map of the Fermi surface. In the former case, 
it maps out the superlattice Brillouin zone generated by 
the nesting instability. 

A number of different mechanisms have been proposed 
for the origin of pseudogaps in the cuprates. These in- 
clude magnetic (spin density waves) ||, flux phase (RVB) 
charge ordering (CDW) ||||, and superconducting 
fluctuations ||. These instabilities fall into two classes 
nesting (the first three: associated with particle-hole 
propagators, and instabilities in the charge or spin sus- 
ceptibilities at a nesting vector, here Q = (%, 7r)) and 



pairing (the last, associated with particle-particle prop- 
agators and the uniform susceptibility at q — 0). In the 
cuprates, all of these instabilities may be analyzed within 
an SO (6) group - the instability group of the Van Hove 
singularity M). We find that the rFs has two strikingly 
different origins in these two classes, but there is rela- 
tively little variation within a given class. In CCOC the 
rFs seems to indicate the locus of the reduced Brillouin 
zone, suggesting a nesting instability. 

An important sum rule for ARPES that relates the 
integrated intensity to the momentum distribution has 
been introduced by Randeria et al. Jilt] and is given by 



n(k) = / duf(cj)A(k,u), 



(1) 



where A(k, uS) is the one particle spectral function of the 
model, n(k) =< ctcj > is the momentum distribution 

and f(ui) is the Fermi function. Although n(k) is a 
ground state property, they proved that in the limit of 
the sudden approximation the frequency integrated spec- 
tral function gives the momentum distribution. They 
employed this sum rule to determine the momentum dis- 
tribution in BSSCO and YBCO. Ronning et al. @ ex- 
tended this methodology to strongly correlated electron 
compounds. By defining kp as the point of steepest de- 
scent, they showed that even when strong Coulomb corre- 
lations destroy the Fermi-liquid character of the system, 
n(k) still drops sharply, allowing the determination of a 
rFs. 

In the following, we calculate the spectral function 
A(k,u>) and the momentum distribution for mean field 
models with a variety of instabilities. Figure [j] illus- 
trates the rFs's associated with several nesting instabil- 
ities. The energy dispersion has the standard one-band 
form 

= — 2fo(cos k x a + cos k y a) — At\ cos k x a cos k y a, (2) 

with to — 0.256^, t\ = —0.45to- The calculations follow 
those in Rcf. |Q for CDW and s-wave superconductivity. 
The rFs's for all the nesting instabilities are essentially 
identical: over part of the surface, there is no gap, and 
the rFs is dispersionless, coinciding with the true Fermi 
surface at = Ep; here n(k) = 1/2 is mainly due to the 
Fermi function in Eq. [I]. Over the rest of the zone, the 
rFs lies along the zone diagonals which determine the 
a/2 x a/2 nesting superlattice. On this jmrt of the rFs 
there is considerable dispersion, and n(k) = 1/2 due to 
the coherence factor, discussed below. By contrast, the 
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rFs for a pairing instability is always located below Ep at 
the superconducting gap, dispersionless for s-wave, dis- 
persive for d-wave, but in both cases faithfully following 
the contours of the true Fermi surface. 
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FIG. 1. Remnant Fermi surfaces for different nesting in- 
stabilities at T—OK. (a) evolution of the rFs toward a perfect 
square with increasing CDW gap; (from darkest to lightest) 
Ocdw=0, 100, 300, 500 meV. (b) the same evolution for a 
flux phase instability. Superconducting instabilities are not 
shown, since in this figure the rFs would coincide with the 
true Fermi surface. 

If the pscudogap is due to a nesting instability compet- 
ing with superconductivity, then there should be a char- 
acteristic evolution of the rFs with doping, from nesting- 
like at half filling to pairing-like in the overdoped regime. 
The phase diagram has been worked out for such a com- 
petition, both for CDW-to-s-wave superconductivity Jl^] 
and for flux phase to d-wave 0]. In both cases, we find 
the evolution of the rFs's is nearly identical. Figure || 
illustrates this evolution for the latter case. Note that 
since the phase at half filling is fully gapped (a Mott in- 
sulator), the rFs is perfectly square. The two limiting 
cases, insulator and optimally doped, bear a marked re- 
semblance to the experimental observations jL|. 

In the calculations of Fig. |[ we took the compet- 
ing phases to be d-wave superconductivity, with gap 
Ai = A d 7£, with 7£ = cos k x a — cos k y a and an or- 



bital antiferromagnot |f3|Jic| ], a nesting instability with 
gap Oi c = jc jj:, which is an RVB state having a d- 
wave symmetry corresponding to a particle-hole excita- 
tion, essentially equivalent to the flux phase instability 
introduced by Affleck and Marston Q]. We consider a 
one-band model, Eq. ||, with correlation effects simulated 
by a doping dependent to = xt$, tg = 2.3eV, and 

tt/to = -0.52tanh(2.4a;), (3) 

to pin the Van Hove singularity (VHS) close to the Fermi 
level over an extended range of doping P]l2| . We start 
with the following mean-field hamiltonian, 

H = erct c,- + Or C (4 ce , a + h.c.) + 

Z-^i « fa " k v ka k+Qa I 

fc,£T 

V A^(ct c f - +h.c). (4) 

k 

with the quasiparticle dispersion given by, 

E h=l(4 + ^ + ^l + Of 2 ±Ei)) (5) 

Performing a generalized Bogoliubov-Valatin transforma- 
tion we find the following gap equations: 

BE 

A rf = U A — ( cos2 sin 2( t ) + tanh — + 

k 

BE - 

sin 2 sin 20_ tanh — ) (6) 

and 

BE - 

O jc = U .jc ^2 -j sin 2( ^ cos tanh — : ~ + 

k 

3E_ jj 

cos 20_ tanh — ) (7) 

with Ej: = J (eg — e k+Q> 2 + ^i° 2 anc ^ tne an gl es de- 
fined by tan20 ± = 2Ag/(eg + eg +<5 ± E % ), 

cos 20= £g ~ £g +Q , ( 8 ) 
E k 

and eg = eg — Ef- Solving the gap equations self- 
consistently with the coupling constants Uo-'c =176 meV 
and Ca=88 meV, we find that the experimental phase 
diagram can be fit semi-quantitatively by this model (in- 
sert in Fig. ||) p5|. This phase diagram differs from that 
found in Ref. Jl2| , but not because different instabilities 
were assumed. The main difference arises from the gap 
cutoffs: in Ref. p2| , a phononic interaction was assumed, 
with cutoff fioj = A5meV; here, the interaction is taken 
as electronic, with no cutoff. This makes a significant 
difference in the phase boundary near half filling. 
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FIG. 2. (a) Evolution of the rFs from half filling to optimal 
doping at T = OiC" with self-consistent gap parameters for a 
model with d-wave superconductivity and flux phase. Curves 
from ligthest to darkest: x=0 (thin line), 0.04, 0.1, 0.19 and 
0.26 (black line). Inset: Pseudogap phase diagram for LSCO 
(filled circles) and YBCO (open squares) determined from 
transport fll5j. Solid line - O jc transition T*; dotted line - 
superconductivity transition T c . (b) quasiparticle dispersion 
along the rFs plotted in (a), tan# = k x /k y 



In the above calculations, the shape of the Fermi sur- 
face evolved with doping, to mimic the effect of Van Hove 
pinning |l6|||. However, the evolution of the rFs for a 
nesting gap is generically toward a square Fermi surface 
as the gap is increased, as illustrated in Fig.l. The rFs 
has two parts, a true (hole pocket) Fermi surface on the 
ungapped part of the rFs, and a segment of square on 
the gapped part. With increasing gap magnitude, the 
former feature shrinks and the latter grows, until the full 
Fermi surface is gapped and the rFs is square. Note that 
due to the d-wave symmetry the rFs for a flux phase in- 
stability is pinned at the true Fermi surface along the 
diagonal (0, 0) — (tt, tt), Fig. lb. The shrinking of the true 
Fermi surface is reminiscent of the evolution in BSCCO 
reported by Norman, et al. fTij] . It should be noted that 
the rFs is not equivalent to the minimum gap locus in- 
troduced by Ding, et al. |Q. 

The origin of the rFs can be understood from these 
calculations. In the competing flux phase-d-wave model, 
n(k) can be written as 



1 (3E g 
n(k) = — (1 — cos cj) cos 20+ tanh — 

sin (f> cos 2(f>- tanh — — — 



(9) 



For a pure d-wave superconductivity model this becomes 
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P E %, 



n(k) = -(l--^tanh ^ 

k 



(10) 



with Et = 



4- 



A|, showing that the rFs coincides 



with the true Fermi surface: n(fc)=l/2 when eg = 0. For 
a pure nesting model n(k) is given by 



n(k) = -(1 



/3E nes J 

> tanh '■ 

2 

(3E n _ es J 
<j) tanh — J — ) , 



(11) 



with E nes J 

±,k 



0, the two 
) — sin 2 (f> = 

e k+Q' ^ or ^ e P resen t model, this 



tanh's go to 1 or -1, so n(k) — 1/2 when cos 2 i 
0, or, from Eq.[8], i 

is the superlattice Brillouin zone boundary. 

In the underdoped regime, as temperature is lowered 
the cuprates pass first into the pseudogap phase, at tem- 
perature T*, then into a superconducting phase at T c . 
In the present scenario, T* would signal a transition to a 
nested phase with a gap (or pseudogap if realistic fluctu- 
ations are included p|,^9[), leaving hole pockets behind. 
Below T c , an additional, pairing gap opens at the hole 
pockets. However, a careful look at the rFs shows a 
more complicated evolution, Fig. 3a: the shape of the hole 
pockets changes, with an accompanying transfer of spec- 
tral weight from the nesting to the pairing parts of the 
rFs. Note that in Fig. 3a the rFs has the same locus in 
k-space as the true hole pocket Fermi surface above Tc, 
but from Fig. 3b there is a dramatic shift in dispersion of 
this rFs as the superconductivity gap opens. 

In comparing these results to experiment, the rFs of 
CCOC clearly displays the square shape expected for 
a predominantly nesting interaction. This is consistent 
with all of the pseudogap models noted above, except for 
preformed pairs. In fact, preformed pairs would still be a 
possibility, if strong correlation effects renormalized the 
(true) Fermi surface to square at half filling. Such renor- 
malization has been proposed previously and is 

incorporated in Eq. |^. However, in these theories, the 
renormalization leads to greatly enhanced nesting, and 
is less favorable for pairing. The most likely conclu- 
sion is that the pseudogap in the underdoped cuprates 
represents some nesting instability, which is fundamen- 
tally competing with superconductivity. Clearly, since 
the cuprates are quasi-two-dimensional, there should be 
prominent superconducting fluctuations above T c , but 
they do not represent the dominant part of the pseu- 
dogap. 
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FIG. 3. (a) Evolution of the rFs with temperature 
for a fixed doping £=0.19: black line - T—0K, grey line - 
T = T C =94K (A d =0 meV), dashed line shows what the flux 
phase rFs at T=0K would be if A d =0 meV. Inset: Tempera- 
ture dependence of the superconducting and flux phase gaps; 
dotted line - A d , solid line - O jc , dashed line - VA d2 + JC2 
(b) quasiparticle dispersion along the rFs plotted in (a). 

The present results suggest a number of experimental 
tests. The rFs should be mapped out in the cuprates as 
a function of doping. In particular, the results of Nor- 
man, et al. |l~7j should be extended to the full rFs. Ob- 
servation of a shift in spectral weight with temperature, 
Fig. |U would provide strong evidence that the pseudo- 
gap is a nesting phenomenon, and not due to preformed 
pairs. Moreover, the rFs can be studied in other systems, 
to confirm the predicted properties. A start has already 
been made in CDW systems p2fl . Mott insulators would 
be particularly of interest. It is believed that the insu- 
lating phase can form in the absence of magnetic order, 
hence without nesting. In the cuprates, there is a clear 
Neel transition, so a square rFs is not unexpected, but 
a study of rFs's in non-magnetic Mott insulators could 
prove most informative. 
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